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Abstract: Accurate automatic identification of astronomical objects in an imperfect world of non- 
linear wide-angle optics, imperfect optics, inaccurately pointed telescopes, and defect-ridden cameras 
is not always a trivial first step. In the past few years, this problem has been exacerbated by the rise of 
digital imaging, providing vast digital streams of astronomical images and data. In the modern age of 
increasing bandwidth, human identifications are many times impracticably slow. In order to perform 
an automatic computer-based analysis of astronomical frames, a quick and accurate identification 
of astronomical objects is required. Such identification must follow a rigorous transformation from 
topocentric celestial coordinates into image coordinates on a CCD frame. This paper presents a fuzzy 
logic based algorithm that estimates needed coordinate transformations in a practical setting. Using 
a training set of reference stars, the algorithm statically builds a fuzzy logic model. At runtime, the 
algorithm uses this model to associate stellar objects visible in the frames to known-catalogued objects, 
and generates files that contain photometry information of objects visible in the frame. Use of this 
algorithm facilitates real-time monitoring of stars and bright transients, allowing identifications and 
alerts to be issued more reliably. The algorithm is being implemented by the Night Sky Live all-sky 
monitoring global network and has shown itself significantly more reliable than the previously used 
non-fuzzy logic algorithm. 

Keywords: data analysis - methods: statistical - techniques: astronomical image processing 



1 Introduction 

Useful pipeline processing of astronomical images de- 
pends on accurate algorithmic decision making. For 
previously identified objects, one of the first steps in 
computer-based analysis of astronomical pictures is an 
association of each object with a known catalog entry. 
This necessary step enables such science as automat- 
ically detected transients and automated photometry 
of stars. Since computing the topocentric coordinates 
of a given known star at a given time is a simple and 
common task, transforming the celestial topocentric 
coordinates to image (x, y) coordinates might provide 
the expected location of any star in the frame. How- 
ever, in many cases slight shifts in the orientation, in- 
accuracy of the optics or imperfections in the CCD can 
make this seemingly simple task formidable. 

Fuzzy logic (IZadeblll965l. Il978l . Il983t 11994) is an 

extension of Boolean logic that is useful for making 
complex decisions by computers. While in Boolean 
logic an item has only two levels of memberships to 
a set (false or true; or 1), fuzzy logic supports any 
value within the interval [0,1] as a level of membership 
of an item to a fuzzy set. A fuzzy logic model consists 
of three different types of entities: fuzzy sets, fuzzy 
variables and fuzzy rules. The membership of a fuzzy 
variable in a fuzzy set is determined by a function that 
produces values within the interval [0,1]. These func- 
tions are called membership functions. Fuzzy variables 
are divided into two groups: antecedent variables that 
contain the input data of the fuzzy logic model, and 



consequent variables that contain the results calculated 
by the fuzzy logic model. 

In this paper, we present a fuzzy logic based algo- 
rithm that transforms celestial coordinates into image 
coordinates for even complex combinations of wide- 
angle non- linear optical distortions, slight optical im- 
perfections, and small unrecorded orientational per- 
turbations. In Section 2 we describe the fuzzy logic 
algorithm, and in Section 3 we present and discuss 
a practical implementation of the algorithm for auto- 
matic analysis of an ope rational astronomical project 
called "Night Sky Live" ijNemiroff et alJl200,^ . 

2 Fuzzy Logic Based Coordi- 
nate Transformations 

The transformation of celestial topocentric spherical 
sky coordinates to local Cartesian image coordinates 
can be defined by a set of two functions. Mathematical 
conversion from right ascension, declination, latitude, 
longitude, and local time to altitude and azimuth is 
straightforward. Further transformation of azimuth 
and altitude sky coordinates to (x, y) image (CCD) 
coordinates is frequently harder in practice when faced 
with practical inaccuracies of knowledge. For practical 
application of fuzzy logic, this later transformation is 
broken up into two parts: 



[altitude, azimuth) 
[altitude, azimuth) 



(1) 
(2) 
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On a CCD image, pixel locations can be specified 
in either Cartesian or polar coordinates. Let x zen be 
the x coordinate (in pixels) of the zenith in the image, 
and y zen be the y coordinate of the zenith, x and y 
coordinates of any given star visible in an astronomical 
image can be computed as follows: 

x — Xz en + sin(angle) ■ distance (3) 

y = t/zcn + cos(angle) ■ distance (4) 

Where angle is a polar azimuthal angle and distance 
is a polar radial distance. 

In order to use these equations it is necessary to 
compute a polar distance and angle for given objects. 
Given the observer's latitude and longitude, this can be 
done by converting the celestial coordinates (azimuth 
and altitude) of a given stellar object at a given time 
to the corresponding angle and distance in the image. 
Since the azimuth and altitude of any given bright star 
or planet at any given time can be easily computed, 
the only missing link here is the transformation of the 
altitude and azimuth to image coordinates, so the ob- 
ject can be found in the image. 

2.1 Reference Stars 

Each of the two models is built based on manually 
identified reference stars. A reference star can be any 
star within an image that was correctly associated with 
the corresponding stellar object. Being familiar with 
the night sky, we can inspect the frame by eye and 
identify the stellar objects that appear in it. The im- 
age (x, y) coordinates of the object can be taken from 
the peak of its point spread function (PSF). Each ref- 
erence star provides a record with the following fields: 
azimuth, altitude, angle, and distance. 

Each identified star contributes an azimuth and al- 
titude (by basic astronomy) and also an angle and dis- 
tance (by measurement from the image). These pro- 
vide the raw data for constructing a mapping between 
the two, using the fuzzy logic model that will be de- 
scribed later in the paper. In order to obtain an accu- 
rate fuzzy logic model that calculates the angle, it is 
necessary to select reference stars that uniformly cover 
the entire image. This assures that the calculation of 
one value will depend on reference points that are rel- 
atively close to it. This is also true for the model that 
calculates the distance. 

2.2 Building the Fuzzy Logic Model 

In order to transform celestial coordinates into image 
coordinates, two different fuzzy logic models are being 
built based upon the two transformations: 

/i : azimuth i — > angle (5) 

/2 : altitude, azimuth i — > distance (6) 

Here, altitude, azimuth and angle are angular mea- 
sures, while distance is measured in pixels. Each trans- 
formation (/i and fz) is computed by a different fuzzy 
logic model, thus one model calculates the angle and 
the other calculates the distance. The fuzzy logic model 
/i, which calculates the angle, is an approximation 



based on the assumption that the optical axis is per- 
fectly aligned, so the angle is not depended on the 
altitude. 

The fuzzy logic model /i has one antecedent (in- 
put) fuzzy variable and one consequent (output) fuzzy 
variable, while the fuzzy logic model has two an- 
tecedent variables (altitude,azimuth) and one conse- 
quent fuzzy variable. 

2.3 Converting Azimuth to Polar An- 
gle on the CCD (f±) 

The first model (/i) is built according to the refer- 
ence stars such that each reference star adds to the 
model one fuzzy set and one fuzzy rule. Each fuzzy 
set is associated with a membe rship functio ns that is 
built in the form of a triangle JZadehlll965l) . Each of 
these membership functions reaches its maximum at 
the reference value, and intersects with the x-axis at 
the reference values of its neighboring reference stars. 
For instance, suppose we would like to build the fuzzy 
logic model with a data set that contains the following 
four reference stars: 
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The first reference star maps azimuth 0°. Assum- 
ing Qi < «2 < olz, the membership functions that will 
be added to the model are described in Figure Q 
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Figure 1: The membership functions of the fuzzy 
sets (FSq to FS4) created by the four ref- 
erence values (0,6 ), (oti,9i), (02,62), (a 3 , 63) for 
fi(azimuth 1 — > angle). 



The membership of the fuzzy sets is determined by 
the following membership functions: 
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F 1 (x) 
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Fo to F4 are the membership functions of the fuzzy 
sets FSo to FS4 that were created according to the ref- 
erence stars (the function F m is the membership func- 
tion of the fuzzy set FS m ). The membership functions 
are built such that almost all azimuth values belong 
(with non-zero membership) to two fuzzy sets. Only 
the points of maximum (0,«i ,0:2,03 ,360) have a non- 
zero membership to just one set. 

The defuzzification method used in this model is weighted 
average, which is an efficient defuzzification method 
when the fuzzy lo gic model is built according to a set o f 
signleton values ST. Takaei fc M. Sugenolll983l Il985h . 
Since weighted average is used, the consequent part of 
each rule is a crisp value and not a fuzzy set. There- 
fore, the fuzzy rules that will be added to the model 
are: 

FSo 1— ► 60 
FSi 1— » 0i 
FS 2 ' — ► 62 

fs 3 1— » e & 

FSi ' — > 360 



For instance, suppose that the first three reference 
stars have azimuths of 0°, 10°, and 20°, and their po- 
lar angles are 3°, 12°, and 22° respectively, such that 
a = 0°, 0o = 3°, Qi = 10°, 0! = 12°, a 2 = 20°, and 
02 = 22°. Each reference star adds one fuzzy set to 
the model such that the membership functions of the 
first two fuzzy sets FSo, FSi are: 



F (x) = 



Fi(x) 



< x < 10 

x < or x > 10 



10 



10) 



< x < 10 
10 < x < 20 
x < or x > 20 



Each reference star also adds one fuzzy rule such that 
the first three fuzzy rules are: 
FSo 1— » 3 
FSi 1 — > 12 
FS 2 1 — ► 22 



Now suppose that we want to use this model in 
order to compute the polar angle (in the image) of a 
stellar object with an azimuth of 6°. The value 6 has 
a level of membership of 1 — • 6 = 0.4 to the fuzzy 
set FSo and ■ 6 = 0.6 to the fuzzy set FSi. Since 
the level of membership of 6 to all other fuzzy sets is 0, 
the only rules that will have any affect on the output 
of the computation are rules and 1 above. Since the 
weighted average defuzzification method is performed, 
the output of the computation is 0,4 3 4 "^° 6 6 12 = 8.4. 
This computation practically provides the same results 
as a simple linear interpolation of the value 6 with the 
two closest neighboring reference stars. 



2.4 Converting Altitude and Azimuth 
to Radial Distance on the CCD 

(/a) 

Unlike the simpler /1 model used for transforming the 
azimuth to angle, the computation of the distance (in 
pixels) from (x zen ,y zen ) should be computed based on 
two parameters, which are the altitude and the az- 
imuth. Using both the altitude and azimuth allows 
the model to deal with asymmetric behavior of the 
optics as well as inaccurate orientational information. 
In other words, using the assumption that the optics 
orientation is directly at the zenith and the distortion 
of the optics and hardware is completely symmetric, a 
reference point at a certain azimuth would allow calcu- 
lating the distance of a point at the same altitude but 
at a different azimuth. However, this is not always the 
case. For instance, a stellar object with the azimuth of 
0° (north) and altitude of 30° can be at distance of 150 
pixels from (x zen ,y zen ), while another stellar object at 
the same altitude (30°) but at azimuth of 60° will be 
at distance of 155 pixels from (x^ en , y zen ) . Moreover, 
(xzen, Uzen) does not necessarily appear in the center 
of the frame, and the frame is not necessarily central- 
ized. Therefore, when computing the distance of a 
stellar object from (a; ze n, Uzen), it is required to take 
not only the altitude of the stellar object into consid- 
erations, but also the azimuth. In order to do that, 
the fuzzy logic model that calculates the distance is 
built according to reference stars that not only have 
different altitudes, but also different azimuths. 

We build the fi model that computes the distance 
using four different sets of reference stars such that 
each set contains reference stars that share approxi- 
mately the same azimuth. For the sake of simplicity, 
the first set contains reference stars that are near az- 
imuth 0°, the second set contains reference stars near 
azimuth 90°, and the other two sets contain stars near 
azimuth of 180° and 270° respectively. I.e., all refer- 
ence stars used for this model should be fairly close 
to the azimuth of 0°, 90°, 180° or 270°. In order to 
use the four sets of reference stars, four new fuzzy 
sets are added to the model. Those fuzzy sets are 
"North", "East", "South", and "West". 

Fuzzy logic can be viewed as a complex interpola- 
tion method. In two dimensions, our fuzzy logic re- 
alization allows /2 to perform the combined interpo- 
lation of linear and Gaussian functions in a straight- 
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forward fashion, and to alter the model easily when 
needed. Note that this is quite different than four sepa- 
rate linear and/or Gaussian interpolations for the four 
directions. The membership functions of the fuzzy sets 
are described in Figure [5] 



Table 1: Reference stars 



North East 
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Vest North 




180 
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Figure 2: The fuzzy sets of the four directions 
(North, South, East, West). 



The membership in the fuzzy sets is determined by 
the following Gaussian functions: 
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90 < x < 360 

< x < 90 

Where a is set to 45. 

In this case, the level of membership for each of 
the four membership functions is always greater than 
zero, as opposed to the /i angle model where only two 
membership functions have a positive level of member- 
ship. 

Building fa model can be demonstrated by an ex- 
ample: Suppose the model is built based on the fol- 
lowing 6 reference stars listed in Table Q 

As with /i, each reference star adds to the model 
one fuzzy set that has a triangle membership function. 
For instance, the membership function of the fuzzy set 
added by Si reaches its maximum of unity at 68, and 
intersects with the x-axis at the points of maximum of 
its neighboring reference stars. The neighboring stars 
are the two stars that their altitudes are closest to the 
altitude of Si (such that one is greater than 68° and 
one is smaller than 68°), and have approximately the 
same azimuth as Si . In this example, one neighboring 
star would be 5*2 and the other would be 5*3 . All three 
stars share approximately the same azimuth (which is 
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North). Therefore, the fuzzy set AU68N that will be 
added by Si will have a triangle membership function 
that reaches its maximum at 68, and intersects with 
the x-axis at 62 and 72. This membership function is 
described in Figure [3] 



Membership 



AlteSN AH72N 




68 
Attitude 



Figure 3: The membership function of the fuzzy 
set AU68N. The figure also include parts of the 
membership functions of AU62N and AU72N that 
were added by the two neighboring reference stars 
S% and S3. 



F, 



The membership function added by Si is: 
~ 62<a;<68 
68 < x < 72 
x < 62 or x > 72 



r(x) = 




Since the azimuth of Si is approximately north, 
the fuzzy rule added to the model by Si is: 
Alt68N A North 1 — > 215 

This rule is more significant for stars that appear in 
the northern part of the sky (and in this case, also at 
an altitude of around 68°). The practical affect of this 
rule will be stronger as the coordinates are closer to 
the 68° parallel. The fuzzy rules that will be added by 
the other reference stars are: 
Alt62N A North 1 — > 224 
Alt72N A North 1 — > 206 
Alt62E A North 1 — > 188 
Alt66E A North 1 — ► 180 
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Alt70E A North i — > 172 

Now suppose that we want to use this model in 
order to compute the distance (in pixels) of a stel- 
lar object with an azimuth of 30° and altitude of 64°. 
The value 64 have a membership value of ggl^ — | 
to the fuzzy set AU68N (that was added to the model 
by ft), and 1 - §§5§§ = § to the fuzzy set AU62N 
(that was added to the model by &)■ Also, this value 
would have a membership value of §§5fl = ^-5 to the 
fuzzy set added by S4, and a membership value of 
1 - §§5§§ = 0.5 to the fuzzy set added by S 5 . The 
membership value of the azimuth 30° to the fuzzy set 

30 2 

North is e TIE 7 ~ 0.8, and to the fuzzy set East it 

(30-90) 2 

is e 2-45 2 ~ 0.412. The membership values to the 
fuzzy sets South and West are very close to in this 
case, and therefore have very little effect on the output 
value. For the sake of the simplicity of the example, 
these membership values are assumed to be exactly 0. 

The computation process is based on product in- 
ferencing and weighted average defuzzification. There- 
fore, the output value of the computation would be: 

215-0. 8-0. 333+224-0. 8-0. 667+180-0. 412-0. 5+188-0. 412-0. 5 ^ 9n o aj 
0.8-0. 333+0. 8-0. 667+0. 412-0. 5+0. 412-0. 5 — 

3 Example Application to Night 
Sky Live Data: Accuracy and 
Complexity 

The fuzzy logic based transformation formula has been 
tested and is i n practical use with t he Night Sky Live 
project (NSL) llNemiroff et al.ll2005fi . Purposes of the 
global Night Sky Live (NSL) network of fisheye CON- 
tinuous CAMeras (CONCAMs) include the ability to 
monitor and archive the entire bright night sky, track 
stellar variability, and search for transients. Fully 2-7T 
steradians - half the sky - are monitored passively by 
each CONCAM, without tracking. Currently, the NSL 
project deploys nine CONtinuous CAMeras (CONCAMs) 
at many of the world's premier observatories. When 
the Moon is down, CONCAMs take 180-second expo- 
sures every 236 seconds, and can detect stars down to 
visual magnitude 6.8 near the image center. 

The fuzzy logic based transformation formula is 
used by NSL for converting the celestial coordinates 
to image coordinates, so known catalogued stellar ob- 
jects can be associated with PSFs that appear in the 
Night Sky Live frames. One simple task that is enabled 
by this transformation formula is the annotation of the 
all-sky images with the names of bright stars, constel- 
lations and planets. This task is mostly used for educa- 
tional or "cosmetic" purposes. A more important task 
is the automatic detection of non-catalogued objects. 
This task is required for automatic detection of me- 
teors, comets, novae and supernovae, as well as other 
astronomical phenomena visible in the night sky. The 
following algorithm uses the transformation formula in 
order to associate PSFs in the image to stellar objects. 

1. function check_stars(image, date_time) 

2. image_PSFs <- GetPSFs(image) 



3. for each image_PSF_cords in image_PSFs do 

4. begin 

5. min_distancc <— 00 

6. for each star in catalog do 

7. star_celestial_cords <— CelestialCoordinates(star, 
date_time) 

8. if InView(star_celestiaLcords) then 

9. image_cords <— F(star_celestiaLcords) 

10. if distance(image_cords,image_PSF_cords) 
< min_distance then 

11. min_distance ^— distance(image_cords, 
image_PSF_cords) 

12. end if 

13. end for 

14. if min_distance < TOLERANCE then 

15. associate(image_PSF_cords,star) 

16. end for 

The function F transforms celestial coordinates to 
image coordinates based on the fuzzy logic models de- 
scribed in this paper. The function GetPSFs returns 
a list of coordinates of the PSF peaks that appear 
in the picture. This function can be implemented 
by using some available algorithms for detection of 
sources from astronomic al images such as SExtractor 
l)Bertin fc Arnouts! ll995L The function InView re- 
turns true if its argument coordinates are inside the rel- 
evant view of the optical device. In the inner loop the 
algorithm searches the catalog for a star that should 
appear closest to the center of the PSF. Since the hard- 
ware used for the Night Sky Live project currently can- 
not get deeper than magnitude 6.8, t he catalog being 
used is a subset of Hipparcos catalog (ESA 1997) that 
is sure to include objects this bright. In line 14, the 
minimum distance found in the inner loop is compared 
with a constant value TOLERANCE that is a tolerance 
value. Only if the distance is smaller than TOLER- 
ANCE then image_cords and image_PSF_cords are 
considered as referring to the same star. 

In Figure g] taken by the Mauna Kea CONCAM, 
the names of bright stars and constellations were an- 
notated automatically using the above transformation. 
The coordinates of one bright point (appears at the 
lower right of the frame) could not be associated with 
any catalogued object so it was automatically marked 
with a yellow square. This object is believed to be a 
meteor. 

The present fuzzy-logic algorithms allow practi- 
cally 100 percent chance of accurate identification for 
NSL stars down to a magnitude of 5.6. We are cur- 
rently unaware of any exceptions. A previously used 
NSL identification algorithm that employed a straight- 
forward analytic transformation was only accurate to 
about magnitude 3.5, although that was somewhat de- 
pendent on the NSL station. This dramatic improve- 
ment was the driving impetus for the creation of this 
paper. 

A useful by-product of the newly accurate iden- 
tifications is the automatic generation of photometry 
files. The ability to associate each PSF with a cat- 
aloged star allows the system to provide continuous 
monitoring of many bright stars. This information is 
provided in the form of text (XML tagged) files. Each 
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frame produces one text file that lists all PSFs that 
were detected in the frame and the name and cata- 
logue number of the stellar object associated with it. 
It also lists some additional data about each detected 
object such as the previously cataloged visual magni- 
tude, spectral type, and celestial coordinates. Identifi- 
cation allows other algorithms to process the frame and 
match each star with real time photometric data such 
as estimated counts of the background and the counts 
of the PSF. For each bright star, the average of the 
brightest 1, 5, 9, 16 and 25 pixels are listed. For dim- 
mer stars, however, the realized PSF is much smaller 
so only the average of the top five pixels is listed. The 
on-line database allows browsing the records for past 
events. Since the photometry data itself is out of the 
scope of t his paper, more infor mation about it can be 
found at jNemiroff et al.ll2005h . 

The inverse computation of the presented transfor- 
mation formula (converting angle and distance to al- 
titude and azimuth) can be built in the same method 
described in this paper, with the exception of using 
the angle and distance for defining the membership 
functions, while the azimuth and altitude are used as 
the crisp output values of the fuzzy rules. This trans- 
formation formula is currently being used by NSL for 
computing 3D trajectories of meteors detected by the 
twin CONCAMs located at Mauna Kea and Haleakala. 

3.1 Computation Accuracy 

In order to check accuracy, the combined /i and ji 
models were used to compute the image coo rdinates of 
some 150 stars recorded by the NSL project IjNemiroff et alJ 
2005). The NSL project uses fisheye images with an 
extreme viewable angle of 180°. For each pair of im- 
age (x,y) coordinates, we calculated the Euclidean dis- 
tance from the location calculated by the model and 
the location where the star appeared in the image. 
We took into considerations the average of the 150 
Euclidean distances (the average error) and the worst 
case error, which is the longest distance among all com- 
puted 150 Euclidean distances. Table [5] below shows 
the accuracy levels according to the number of refer- 
ence stars that were used in order to build both fuzzy 
logic models (the altitude /2 model and the azimuth 
/i model). 



Table 2: Accuracy Level of the Transformation as 
a Factor of the Number of Reference Stars 
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3.2 Computation Complexity 

The complexity of the computation is a function of 
the number of fuzzy rules and fuzzy sets in the model, 
which, in turn, is depended on the number of reference 
stars. The accuracy needed to identify stars unambigu- 
ously is depended on the density of stars per pixel in 
the frame. Clearly, the higher the average star density, 
the greater the accuracy needed to avoid confusion. 

Experiments suggest that accuracy within five pix- 
els will allow automatic analysis of Night Sky Live 
frames. As shown in Table |21 this can be achieved 
using a fuzzy logic model that is built with constant 
number of around 80 well spread out reference stars. 
Therefore, the final fuzzy logic model contains a con- 
stant number of around 80 rules and around 160 fuzzy 
sets, so the theoretical complexity of the computation 
is O(l). Practically, the CPU time that is required for 
the computation is negligible, and around 1000 coor- 
dinates are transformed in less than one second by a 
system equipped with an Intel Pentium III processor 
at 850 MHZ and 128MB of RAM. 

4 Conclusions 

Inaccuracies of optics and hardware lead to non-trivial 
transformation formulae that are sometimes required 
for the automatic analysis of the digital frames. A 
fuzzy logic based method does not achieve an exact 
solution to the problem, but rather provides simple 
and maintainable transformation formulae with an ac- 
curacy that makes them more useful in practice than 
previously attempted straightforward analytic trans- 
formations (based on the linear transformation R = fee, 
where R is the radial distance from x zen ,y zen , e is the 
altitude and k is some constant). When applied to the 
NSL project, the presented method is accurate enough 
to perform automatic processing of all-sky images in 
order to associate PSFs in the frames with their cor- 
responding stellar objects nearly 100 percent of the 
time down to magnitude 5.6, and automatically detect 
non-catalogued bright objects. This technique demon- 
strates that it is possible to use fuzzy logic based al- 
gorithms to reduce the affect of minor defects and in- 
accuracies of the optical hardware and/or slightly in- 
accurate telescope orientation information. For NSL 
frames, only 80 reference stars are required to build 
the fuzzy logic model so that automatic identifications 
do not form a computational bottleneck. 
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Figure 4: A Night Sky Live picture with labeled stars, constellations, a planet and a meteor (boxed) 
processed using the fuzzy logic based transformation formula. 



